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Abstract 

The method of the statistical sample moments was used for the analysis of 
neutrino events from SN1987A burst in Cherenkov detectors. In particular the co- 
efficients of correlation Q(E, t) between the energies E of electron antineutrinos 
v e emitted by star and the ejection instants t of v e for neutrino events recorded 
by Cherenkov water detectors of KII and 1MB collaborations were calculated. 
Q(E, t) values depend on the assumed mass of v e . Modern model of the grav- 
itational stellar core collapse with an accretion phase predicts the low level of such 
correlation (Q(E,t)) averaged with respect to neutrino burst. On condition that 
empirical Q(E, t) values equal the theoretical model quantities {Q(E, t)) one can 
obtain 22 ± 10 eV/ c 2 as an estimate of the nonzero v e mass. The error of this 
estimate implies that v e mass less than 2eV/c 2 is unlikely. The laboratory data 
of the tritium /3-decay agree adequately with the presented astrophysical estimate 
provided that the anomalous structure near the end point of /3-spectrum is taken 
into account. 

PACS: 14.60.Pq, 97.60.Bw 

1 Introduction 

The observation of SN1987A supernova neutrino burst gave valuable information re- 
lated to the star collapse astrophysics and neutrino physics. Mainly this information 
was obtained by Cherenkov water detectors of Kamiokande II (KII) [1] and Irvine- 
Michigan-Brookhaven (1MB) [2] collaborations. Two decades before the burst, Zat- 
sepin suggested the possibility of estimation of the neutrino mass taking into account 
that massive neutrino with rest energy E and total energy E passes distance d from 
the source to detector in a time [3] 
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Here c is the velocity of light. Large number of works was dedicated to the esti- 
mation of Eq beginning from 1987, in which as a rule, the parameters of individual 
neutrino events, namely E and the time of neutrino registration ta were compared. 
Such analysis has large uncertainty due to the finite duration of neutrino burst. An 
opinion was formed that the results of observations agree with zero neutrino mass or 
allow one to estimate only its upper limit. 



2 Statistical analysis 

In this paper the results of estimation the electron antineutrino v e mass by comparison 
of the lower empirical sample central moments of two dimensional density function 
F(E, t) of the source and of similar moments of this function predicted by a modern 
theoretical model [4] are presented [j] Empirical moments of F(E,t) depend on the 
assumed value of Eq because td = t + t s d, where t is the instant of the D e emission 
from the source. As usual, the energy of neutrino during its motion from the source to 
detector is assumed to be constant. 

Density function at distance d from the source is 

/>00 

F(E,t d )= / S(t d -t-t sd )F(E,t)dt (2) 
Jo 

Here zero time (t = 0) corresponds to the beginning of emission of v e . If one con- 
siders a specific detector, it is convenient to normalize F(E, t) in formula Q to unity 
in domain < t < oo and E c < E < oo. Here E c designates energy threshold of v e 
detection in considered detector. All sample empirical moments of F(E, t) were cal- 
culated with weight coefficients a in order that they can be compared with theoretical 
ones. For i-th event in the sample of n events 

n 

a, = Eff- 1 (E + )(n- 1 Y / E ^fr 1 (E + ))- 1 (3) 

i=l 

Such structure of corresponds to the process of inverse /3-decay, 

v e + p — > n + e+ (4) 

which, as it is assumed, causes the registration of v e in Cherenkov detectors. Really, 
process Q is 15-20 times more probable than the competing process of the scattering 
of neutrinos (of all types) on electrons. Besides, precisely Q allows one to com- 
pare the energy of neutrino and the number of events observed in [1,2] with predic- 
tions of the stellar core collapse theoretical models. Formula (|3]l takes into account 

1 Below F{x,y) and Q(x, y) stand for respectively two dimensional density function and the coefficient 
of correlation of variables x and y. Furthermore, D(z) and S(z) are the variance and the standard deviation 
of quantity z respectively. 
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the quadratic growth of the capture cross section Q with energy E and registration 
probability f(E+) of a positron with energy E + in a detector (data of KII and 1MB 
collaborations presented in [5] were used in computations). 

Usually the relation between E and E + is expressed by formula 

E = E + + l.3MeV (5) 

In this work a more accurate relation [6], which takes into account a small cor- 
rection related to the recoil energy of neutron is used. However, in the case when 
E + > E c , expression |3]l is sufficiently accurate and it is not necessary to use informa- 
tion about the angular distribution of positron tracks in detectors for process If one 
considers elastic scattering of neutrinos on electrons it is necessary to use such (very 
inaccurate) information. 

The properties of the events from neutrino burst in KII [1] and 1MB [2] detectors 
are summarized in the Tab. [2] 

Table 1 : Measured properties of the positron events detected KII and 1MB in the neu- 
trino burst" 





KII detector [1] 




1MB detector [2] 


Event 


Event 


Positron energy 


Event 


Event 


Positron energy 


num- 


time td 


E+ (MeV) 


number 


time td 


E+ (MeV) 


ber 6 


(sec) 






(sec) 




1 





20.0 ± 2.9 


1 





38 ± 9.5 


2 


0.107 


13.5 ±3.2 


2 


0.42 


37 ± 9.3 


3 


0.303 


7.5 ± 2.0 


3 


0.65 


40 ± 10 


4 


0.324 


9.2 ± 2.7 


4 


1.15 


35 ± 8.8 


5 


0.507 


12.8 ± 2.9 


5 


1.57 


29 ± 7.3 


7 


1.541 


35.4 ± 8.0 


6 


2.69 


37 ± 9.3 


8 


1.728 


21.0 ±4.2 


7 


5.01 


20 ±5 


9 


1.915 


19.8 ± 3.2 


8 


5.59 


24 ±6 


10 


9.219 


8.6 ±2.7 








11 


10.433 


13.0 ±2.6 








12 


12.439 


8.9 ± 1.9 









"It is suggested that all events were caused by the capture of i/ e . 
6 Event no.6 has been excluded by authors from the signal analysis. 



The distance to SN1987A is probably 55 ± 15 Kpc [6]. The value d = 56 Kpc was 
chosen as well as in [4]. When we carry out the statistical analysis of neutrino events, 
it is significant that they all should be caused by the capture of D e (j4j). According to the 
estimates of authors of papers [1,2], the background count rate during 10 s was ~ 0.2 
(KII) and ~ 0.8 (1MB). One event in 1MB was eliminated by authors as a background 
one. I considered all other 8 events. Amongst those eleven events presented by KII 
event no. 3 has value E + = 7.5 ± 2.0 MeV, which practically coincide with energy 
threshold of KII detector. With large probability this event may be considered as a 
background one [6], and it was not taken into account in our computations. Thus, two 
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statistical samples of events (E,td) were taken as initial data, one with size n = 10 
(KII) and another with n = 8 (1MB). Each sample has its own function F(E, t; E c ) 
as population, since values of E c (as parameters) for considered detectors differ es- 
sentially. According to for each of the assumed values of E , the values tj were 
determined from values recorded for considered detector, which gave the sample of 
n events (-E 1 ,, U) specifying density function F(E, t; E c ) of the source. Sample vari- 
ance D(t) and correlation coefficients Q(E,t) and Q(E~~ 2 ,t) were computed. Fig. |T| 
and Fig.|2]represent the dependencies of these quantities on Eq, 

Let us note that these sample moments are invariant with respect to time shift. It is 
essential for our problem and makes possible comparison of the results obtained from 
detectors which are not synchronized in time. Let us emphasize also that considered 
functions do not depend on astrophysical model (except the assumption that process of 
capture Q is real). It is evident from figures that curves D(t) have minima. 

Is it possible to associate the location of the minimum on Eq scale with sought 
value of the v e mass? Let us consider the relation between D(t) on the one hand and 
Q(E~ 2 , t) and Q{E, t) on the other hand. It can be demonstrated from ^ that 

D(t) = D(t d ) - |<f 2 e" 2 • D(E- 2 )E* - dc- 1 ■ Q(E~ 2 ,t) ■ S(E- 2 ) ■ S(t) ■ E 2 (6) 

where S(t) and S(E~ 2 ) are the standard deviations of quantities t and E~ 2 . The first 
term in (|6ji is the variance of the observed values of t c i (it coincides with D (t) when 
Eq = 0), the second term is proportional to the variance of E~ 2 values, and Q(E~ 2 , t) 
is the coefficient of correlation of t and E~ 2 variables in the third term. If one uses 
sample values of statistical moments in ([6}, relation |6) produces just the same values 
of D(t) as a direct method presented above. Expression (|6| allows to understand the 
behaviour of curve D{t) in dependence on Eq. When values of Eq are small, two last 
terms are negative (Q(E~ 2 , t) > 0). When Eq increases, the last term becomes pos- 
itive, since inequality (Q(E~ 2 ,t) < 0) becomes correct. Computations demonstrate 
that the minimum in considered curve precisely corresponds to such value of Eq where 
Q(E~ 2 , t) changes its sign. In this region Eq both functions Q(E~ 2 , t) and Q(E, t) 
have large Eq derivative. The aforesaid is illustrated by Fig. [T] and Fig. [2] It is also 
evident from figures that standard correlation coefficient Q(E,t) changes its sign at 
approximately the same value of Eq as Q(E~ 2 , t). Just function Q(E, t) is considered 
below. It is evident from Fig. [T] and Fig. [2] that functions Q(E,td) corresponding to 
direct observations either are close to unity with respect to modulus (1MB) or have the 
same order of magnitude (KII), i.e. there evident correlation "E — t^" exist. Func- 
tion Q(E, td) specifies "E-time" correlation near detector. Only modern astrophysical 
models of stellar core collapse can present information concerning Q(E, t) near the 
source. 

The results of [4] were used in order to compute coefficient of correlation (Q(E, t)) 
averaged with respect to neutrino burst. In [4] the equations of hydrodynamics (in the 
framework of the general theory of relativity) and equations of multigroup (with re- 
spect to energy) transport of neutrinos which takes into account all types of neutrinos 
and antineutrinos were solved jointly numerically. In particular in [4] luminosity for 
v e and mean instantaneous (at t = const) energy of v e are represented as functions of 
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time which are necessary for computation of (Q(E, t)). At the stage of the core cooling 
(t > 1.2s) exponential decrease of the neutrinosphere temperature was assumed with 
time constant t t = 5.6s in accordance with graphical data [4]. Momentary energy 
spectra of v e were approximated by Fermi-Dirac (F-D) distributions. Such spectra 
which were computed in [4] differ from F-D spectra. They demonstrate steeper de- 
crease at high energies of T> e . The model energy spectra of v e with the same peculiarity 
at high energies [7] were used in one version of our computations. In this case resulting 
values of (Q(E, t)) differ less than to 10%. Possibly it is caused by the fact that mean 
momentary energies of D e were fixed in computations according to data of [4] . The fol- 
lowing approach was used in order to compute the model theoretical value of (Q(E, t)) 
in the same manner as in the case of the sample empirical values of Q(E, t). Values 
of Ei and ti were computed by Monte Carlo method in order to obtain samples of size 
n (as a rule, these sizes coincided with empirical ones). Further value of (Q(E,t)} 
was computed for each k-th sample, and averaging with respect to k(k max = 500) 
was carried out. Obtained strips of values of (Q(E, t)) limited by point straight lines 
as a result of the variation of energy spectra of v e and sizes n are presented in Fig. [T] 
and Fig. [5] It is evident from figures that modulus of (Q(E, t)) is much less that unity, 
i.e. the level of "E — t" correlations in neutrino burst is small on the average. This 
is caused by a large role of accretion phase in the process of gravitational collapse of 
stellar core. For cores with masses k, 2Mq more than half of v e are emitted at the stage 
of accretion specified by positive sign of Q(E, t) [4,5]. Phase of cooling corresponds 
to negative values of Q(E, t). Thus both phases to some extent compensate each other 
in computation of (Q(E, t)). High energy threshold E c of neutrino detection increases 
the role of accretion phase since v e with higher energies are emitted at this phase. Val- 
ues E c = 20 MeV (1MB) and E c = 8 MeV (KII) were used in computations. Thus 
one can understand a positive gap between the values of (Q(E,t)) for 1MB on the 
one hand and of KII on the other hand (see Fig. [T]and Fig. [2]). Large (with respect 
to modulus) and negative (with respect to the sign) differences (Q(E, t) — (Q(E, t))) 
specifying neutrino events in two detectors, evidence that observed coefficients of cor- 
relation Q(E, td) are mainly related to the effect of the delay of the neutrinos with less 
energy according to ([T]), i.e. suggest nonzero mass of v e . 

Determining the curves Q(E, t) intersection regions (as functions of Eq) with strips 
(Q(E, t)) we arrive at the estimates of E as 29 -=- 31 eV for 1MB and 14 16 eV for 
KII. If we assume that the accuracy of 1MB and KII measurements are equal, Eq of 
electron antineutrino can be estimated as 

E o = 22±l0eV (7) 

Here presented error is the root-mean-square deviation for the sample of the size 
n = 2. Assuming that this error has random nature, let us employ Student's distribution 
with degrees of freedom m = 1 for the estimation of the probability of the assumption 
that E < 2 eV is valid. This probability does not exceed 0.15. Thus we can discard 
this hypothesis at 85% significance level. 

If we want to answer the question formulated above, one may say that the minimum 
of D(t) (according to Eq scale) corresponds to the region of the possible values of Eq 
since this minimum is related to the low level of "E — t" correlations. 
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Aside from random errors, "true" value of Eq can be distorted by initial assumption 
that all neutrino events are caused by process of capture Q. In [6] an assumption was 
considered that event no. 4 in KII was caused by the scattering of neutrino on electron. 
If one takes account of such possibility, it complicates the analysis since in particular 
it requires the construction of function F(E, t; E c ) including not only v e , but all other 
types of neutrinos and antineutrinos. Nevertheless let us consider for illustration how 
the values of D(t) and Q(E, t) vary (depending on Eq) if we interpret event no.4 (KII) 
as the scattering of v e . Considering gravitational collapses of stellar cores one limits the 
energy of emitted neutrinos by ~ 50 MeV . In this case the recoil electrons with energy 
exceeding E c should be scattered through angles ip < 20°. If to assume that cp = 10° 
for event no.4, we arrive at E = 29 MeV (instead of 10.6 MeV in the case of capture). 
Such value of Ei(i = 4) leads to essential variation of the D(t) and Q(E, t) values. In 
particular, D(td) — 24.6(27), minimum value of the variance D(t m i n ) — 9(20) and 
Q(E, td) = — 0.68(— 0.42). The numbers in parentheses correspond to the process of 
the capture. The intersection of Eq axis with curve Q(E, t) takes place at the value of 
Eq which is larger by 3 eV than the value represented in Fig. [2] For the estimation of 
Eq in this case we obtain 

Eq = 24 ± 8 eV (8) 

Considered example demonstrates that taking into account the possibility v e scat- 
tering, allows one to make more similar the results obtained by data from both detec- 
tors, however it does not change essentially the estimate of P e mass. 

We used errors 6E (standard deviations) [1,2] in order to take into account the 
influence of inaccuracies in the determination of positrons E + energy in Cherenkov 
detectors. We computed the values E +i for Gaussian distribution with prescribed value 
of 5E + i and M E + i = E + i by Monte Carlo method. The results of computations were 
the samples with sizes n ( (n = 10) (KII) and n = 8 (1MB)) for E + and E. The 
coefficient of correlation Q(E, t) for a number of values of Eq was computed in each 
sample of values of E. A number of k samples (k = 500) were used for a fixed value of 
Eq in order to obtain a distribution (histogram) according to the value of Q. Examples 
of such distributions are presented in Fig. [3] 

Dashed lines in Fig. [T] and Fig. [2] limit FWHM of these distributions (smoothed 
histograms). It is evident that the distributions have sufficiently large FWHM. The 
point is: how significant is the difference of the values of Q(E,t) for Eq = -j- 
2 eV (the carried out analysis cannot distinguish Eq in this interval) and for the value 
of E — 22 eV estimated the mass of T> e by the method presented above? Fig. [3] 
represents distributions with respect to Q for Eq = 22 eV and similar distributions 
for Eq = 0-^2 eV. The regions of overlapping allow one to estimate the probability 
W of coincidence of the distributions for Eq = -j- 2 eV and Eq = 22 eV in the 
following manner: W(IMB) = 0.54 and W(KII) = 0.07. The probability of the 
fact that despite the data of both detectors (independent measurements), distributions 
for E = + 2 eV and Eq = 22 eV do not differ is 

W = W(IMB) x W{KII) = 0.04. 
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Thus the statement that our analysis cannot distinguish Eq = -f- 2 eV and the 
estimate of Eq in (|7]l can be rejected at significance level not less than 90%. 

In the analysis presented above we used data [4] related to an iron stellar core with 
mass 2M & . Smooth curve in Fig. |4] represents integral neutrino event rate for 1MB 
detector which was theoretically predicted [4]. Step integral spectrum represented by 
dotted line corresponds to experimental data. It is evident that these spectra do not 
agree well. On this basis author of [4] concluded that one should pass to the core of 
different mass (1.35M Q ). Indeed, with the employment of Kolmogorov-Smirnov cri- 
terion one can reject theoretical model for 2Mq at significance level ~ 99%. However, 
the time delay of slow neutrinos with respect to the fast ones on the path from the 
source to the detector is not taken into account in [4], i.e. the case Eq = was im- 
plicitly considered. Solid step lines in the same figure represent integral spectrum for 
Eq — 30 eV. This spectrum much better agrees with theory and one can demonstrate 
that in this case theoretical model of the core with mass 2M® cannot be rejected with 
usually used significance levels. Let us note that value Eq = 30 eV determined by 
fitting from the condition of best agreement with theoretical curve coincide with value 
obtained by the method presented above. 

It is interesting to compare estimates (|7]i and ([8]) with experimental results where 
the mass of v e was determined from data of tritium /3-decay. New result of these 
experiments is the discovery of anomalous structure of /3-spectra [8, 9] - a very weak 
line near end point of /3-spectrum. In [10, 11] an approach was proposed which takes 
into account the anomalous structure and leads to the value Eq = 20 ± 5 eV. 

This result agrees with estimates (|7]i and p). 

3 Conclusions 

Modern model of the gravitational stellar core collapse with an accretion phase predicts 
the low level of correlation between the energies E of electron antineutrinos v e emitted 
by the star and the ejection instants t of v e on the average. 

Similar coefficients Q(E, t) of "E — t" correlation can be calculated for neutrino 
events recorded by Cherenkov water detectors of KII and 1MB collaborations from 
SN1987A neutrino burst. 

These values Q(E, t) depend on the assumed mass of V e . On condition that empiri- 
cal values Q(E, t) must be equal to the theoretical model quantities (Q(E, t)) averaged 
with respect to neutrino burst it can obtained the value ~ 20 eV/c 2 as an estimate of 
the nonzero v e mass. 

The error of this estimate allows to consider that the v e mass less than 2 eV/ c 2 can 
be realized only with little likelihood. 

There are two ways of the v e mass estimation: one from the astrophysical observa- 
tions by the method presented above, and the other from laboratory data of the tritium 
/3-decay. Both methods agree adequately when taken into account the anomalous struc- 
ture near end point of the /3-spectrum. 
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Figure 1: The results of the statistical analysis of neutrino events in 1MB detector, 
(a) The dependence of variance D(t) on E (assumed rest energy of D e ). (b) The 
dependence of the coefficients of correlation on E n : Q(E~ 2 , t) - curve 1 and Q(E, t) 
- curve 2. Dashed lines limit FWHM of distributions according to Q(E,t) obtained 
by Monte Carlo method. The straight dotted lines designate the strip of the model 
theoretical values of (Q(E, t)) averaged with respect to the neutrino burst. 
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Figure 2: The same as in Fig.[T]for neutrino events in KII detector. 
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Figure 3: Distributions of Q(E, t) computed by Monte Carlo method for detectors 1MB 
(a) and KII (b). The errors of determination of the energy of positrons in detectors are 
taken into account. Thin lines represent distribution for Eq = -j- 2 eV, thick lines - 
for E = 22 eV. The regions of the overlapping of histograms are dashed. 
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Figure 4: Integral rate of neutrinos for detector 1MB: smooth curve corresponds to 
theoretical model [4], step spectra correspond to real neutrino events (dotted line cor- 
responds to E — eV, solid line - to E = 30 eV). 
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